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ABSTRACT 

The article presents an analysis of the deformable shell behavior of a road roller, the compacted material under the 
compacting roller of the road roller, in which the rigid circular shell of the roller is replaced by a forcefully deformable 
elliptical shape, which, unlike the circular design, allows variation, adjust and optimize the impact of the road roller on 
the material to be compacted. It is necessary to clarify the general conceptual approach adopted in the first part of this 
article by a nonlinear mechanical-mathematical model of elastic gently sloping arched system under the conditions of its 
cylindrical bend using the curvilinear coordinate with a confirmed simplification about the constancy of the 

curvature of the curve. 
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1. INTRODUCTION 

The mathematical interpretation of the stress-strain state of a cylindrical shell element using differential equations, 
considered in previous studies [1-4, 21], is an approximate one, since it is based on a flat-plate model having an 
initial death [4]. This simplifying assumption did not allow, with sufficient objectivity and reliability, to quantify 
the possible stability loss of the shell in the local area of its contact with the road surface material to be compacted. 
This assertion-assumption is legitimate from the point of view of the almost complete analogy of the deformation 
shell of a double-hinged arch under the action of a uniformly distributed load q = where, even in the linear 

formulation of the problem, the upper level of the critical pressure 17 ^ 1 ;is determined [4, 21]. 

Therefore, it is necessary to clarify the general conceptual approach adopted in the first part of this article 
by a nonlinear mechanical-mathematical model of an elastic gently sloping arched system under the conditions of 
its cylindrical bend using the curvilinear coordinate with a confirmed simplification about the constancy of the 

curvature of the curve, and also provided that the external loads ijp = and deflection m = J, there are 

directed along the normal to the initial surface yq (r [4, 21] (Figure 1). 
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2. MATHEMATICAL MODEL 

Considering the nature of the pressure function and the fixing scheme of the element shown in Figure 1 should be 

considered that the subcritical compressive stress corresponds to the main moment state of the shell (Figure 2), in 

which it is necessary to take into account the bending forces [4]. For this purpose, we introduce, as in the linear theory 
instead of the force parameter i^p, a fictitious transverse radial load [4, 21]. 






(1) 


equivalent to the action of stress and equal to the projection of the distributed force aox to the direction of the 


normal to the curved surface (Figure 1). From equilibrium condition [5-9, 21] 


_. 1 _r ^ ~ 

' dft J ' & h-6 ij ■“ ” ' djft . 
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( 2 ) 



Figure 1: Refined Theoretical Design Model of the Shell 
Element in the Form of a Flat Cylindrical Shell-Arch [4] to 
Assess its Local Stability. 


Replacing iJb = "r —= according to (1), (2), we obtain a refined-transformed nonlinear 

T V ' y fp 


differential equation 

dx\ V-djirf irJ3.r / ^ >J Tp 2 -e(^']>v dx\ 2-E(^)..’j.-f^ 


where the constant A is found using the same method, but using the relative deformation formula 
with curvilinear coordinate [4, 21] when counting the movement parallel to the arc (Figure 1), that is 
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(4) 


after substituting co and performing the integration procedure, taking into account compliance with the kinematic 
boundary condition = 0: 


= «Ul3 = j 


1 -(—X -r 

E 2 ttR^ 


A , fxl xj \ 




Revealing obvious boundary equalities uC03 = 0 and it(il|j) = 0, it is obtained (figure 1): 


(5) 


J.5. \7 irJ3^ C/ 


( 6 ) 


The equation (2) is calculated using the Bubnov-Galerkin method [4, 5, 9, 10], using the already known values (6) 
and tables of type integrals [10-12, 21]: 

2 CxJ . (l - ^ rij:, = 0. = (7) 

f'fi= £4/ /Ifty .4f \ J Sf 257r'/ 

12(1 - S Efs “ llS J ?Eb ■ wR, 0 ? E_b ■ tiR, r “ 64 ■ J ' 2FCfl ^ 

. _Ei?_ _n- (Q\ 

"^i:2 Q-u^3' rj ' u>5rp ' 2 -EB;f 4-T^ ‘ 

whereV^(rj^) - the left side of expression (2) in the form of a differential-algebraic operator and an additional 
definite integral [12, 21] 




££i+£il fir 



= 4f 



7 15 \ . r, 
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64 Ip ■ 


(9) 


After the term-by-term division of expression (8) into / and simple transformations, we have, in a finite explicit 
form, the characteristic second-order functional dependence [13-14] between q and/: 
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+ /r3pyil] + + 


( 10 ) 


At / = 0 we obtain from equality (10) the general design formula for the upper critical pressure [4] 

which is for round (f = 0. = 1,5708 [2]] and elliptical ({ ={_£, = 0,57, ffCO^S?) = 1,434 [2]] shells is 

written as follows: 


(K)_LffS4 £ [ ^ j_MVM 

a75u: ^ R. 3 WsJ J3|j" 

[E) _ SM® „ £ f ^ ^ /L&\ n -1 Anss -p’ 

a75ir 13. Usi ■' J3|. 


0..37S04. 


Using the derived relations (12), (13), determine 

I™ = 50 mm, R. = 600 mm, i5 = S..S mm., E = 156000 f—,y- = 0,2S6= [2-7], tt = 3.1416: 
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Vj.57 J 


R-f ' r r ' -21 + +(-)' & ■ (—)4 ■ 0,37504 = 12,4 


The stability condition is observed with a large margin, as in accordance with the calculated data (14), (15). 

= iS=«»k] = L5. 


(14) 


(15) 


4""^ = & = ^=» K] = 




(16) 


Equating to zero the first derivative ^ = t) function (104), determine the maximum deflection ^ (Figure 1) 


j _ n ^ _ 

“ Wb/ it vt - ^ -(J__ 


IPS ^ 


!§, £(f5 

adequate lower practical load^jj^ [4], including (at tt = 3,1416): 

- for shell of circular profile (^ =0, E’en) = 1,5708) 

„ - «! , t _ 1^1 , -hr™™ 

\L6/ -6W Z LZE [J.-(®,Z5*,''2]W^ 

- for elliptical shell (f = = 0,57, FC0,S7) = 1,434 [2]) 


(17) 


(18) 
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AE) _ (^\ _ LD5 _ U,5.1^^ , 14Lfe _ 


(19) 


The derived formulas (9) - (12), (17) are illustrated in Figure 2 of the existing classical curves [4] (see Table 1) 


Table 1: Values of the Function (20) and (21) to the Construction of the 
Corresponding Graphs presented in Figure 2 


/,mm 

0 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 


21,37 

16,43 

12,74 

10,32 

9,2 

9,31 

10,7 

13,37 

17,38 

22,52 

29 

J!L 

12,4 

9,74 

8,04 

7,29 

7,49 

8,64 

10,74 

13,8 

17,8 

22,76 

28,67 


jjW ^ = 3691,186(0,000043/= - 0,0007568/ + 0,00579), (20) 


^ ^ 3691,186(0,0000322/= - 0,0004235/ 4 0,003356), (21) 


by analogy with the cubic dependencies of the simplified model. The graphs of the functional relations (20), (21), 
vividly characterizing the equilibrium state of the refined-modified mechanical system (Figures 1,2), along with the regular 
values of table 1 contain special points including upper and bottom extreme critical pressures 


Pf iE> _ H T , iv 
7 ■■ f/tlT — 14.-47 J J . 


= 21,37 

mm.' 

9e[p “ 9^00 2 ^ 7 


as well as corresponding movements 

= /*‘^ = 0. ] 

= 8, 8mh, = 6,58 mm. J 


( 22 ) 


(23) 


It is also marked by dotted thickened lines and the maximum possible operational extremes. 


^ mm * 






1 

X = 6.57 1 

12 

7 6 

f. ////// 


13.16 



1 - Function= q^'^(f]for round shell (20) 

2 - Functionq = q ^-^'^(fyor elliptical shell (21) 

Figure 2: Diagrams of "Load-Deflection” Corresponding to 
the Modified-Modified Functional Relationships (20) and (21) 
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Table 2: Information on Loads, Material Parameters (Steel), 
and the Results of the Model Breakdown into Cells 


No. 

Indicator 

Value 

JVo 

Indicator 

Value 

1 

Yieldpoint [MPa] 

235 

10 

Pressure: 1 

700 N 

2 

Elasticmodulusnormal [MPa] 

200000 

11 

Pressure: 2 

700 N 

3 

Poisson'sratio 

0.3 

12 

Pressure: 3 

200 N 

4 

Density [kg / m3] 

7800 

13 

Itemtype 

4nodaltetrahedra 

5 

Temperature coefficient of linear 
expansion [1 / C] 

0.000012 

14 

Maximum length of the side of the 
element [mm] 

4 

6 

Thermal Conductivity [W / (m * 
C)] 

55 

15 

Maximumsurfacecondensationratio 

1 

7 

CompressiveStrength [MPa] 

410 

16 

Dilutionratioinvolume 

1,5 

8 

EnduranceLimit [MPa] 

209 

17 

The number of finite elements 

195019 

9 

Torsionendurancelimit [MPa] 

139 

18 

Numberofnodes 

64262 


The boundary conditions are specified with regard to symmetry, according to which the points located in the cross 
section of the cylindrical shell, located at a distance of half the length of the cylinder from any of the ends, cannot move 
along the z axis. The points located on the end of the shell are prohibited to move in the Oxy plane. Points located on a 
circle lying in the plane of symmetry are not allowed to move along the z axis [21]. 

Further, to solve the problem, a grid with finite elements was built (Figure 4). The results are contained in tables 2 
and 3, and are reflected in figures 9-16. 




a b 

Figure 3: Experimental Model of a Roller of a Road Roller with an Elastic Shell, 
Capable of Forcibly Changing its Shape: a) Partial Assembly without 
Unclamping Rollers; b) Complete Assembly. 



a b 

Figure 4: Model of an Elastic, Thin-Walled Drum (a). Divided into Cells by a Finite-Element Mesh (b) 
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Table 3: Inertial Characteristics of the Model 


No. 

Indicator 

Value 

1 

Mass of model [kg] 

14.518466 

2 

Center of gravity of the model [m] 

(-0 ; -0 ; -0.000002 ) 

3 

Moments of inertia of the model relative to 
the center of mass [kg * m2] 

( 0.566375 ; 0.392075 ; 0.392075 ) 

4 

Reactive moment relative to the center of 
mass [N * m] 

(-0.000049 ; 0.000045 ; 0 ) 

5 

The total reaction of the supports [N] 

( 0 ; 0 ; -197.880822 ) 

6 

Absolutereactionvalue [H] 

197.880822 

7 

Absolute moment value [N * m] 

0.000066 


Consider the results of static machine calculation in the figures. 


11.21 
10.53 
9.846 
9.164 
8 482 
7.8 
7.118 
6,436 
5.754 
5.072 
4.39 
3.708 
3.026 
2.344 
1.662 
0.9801 
0.2982 

Figure 5: Distribution of von Mises Equivalent Stresses. 
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0 2822 
0268 
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0.1482 

0.1276 
0.1135 
0.09927 
0.08509 
0 07091 
0.05673 
0.04254 
0.02836 
0.01418 
0 

Figure 6: Total Linear Displacement. 
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Figure 7: Action on Stress Safety Factor. 
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Table 2: Information on Loads, Material Parameters (Steel), 
and the Results of the Model Breakdown into Cells 
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circle lying in the plane of symmetry are not allowed to move along the z axis [21]. 
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Figure 11: 5th form of Natural Oscillations of the Shell 

The critical efforts obtained theoretically are several times larger than the values obtained using modeling in the 

APM FEM KOMPAS-3D, and therefore (as it happened) the results obtained using numerical analysis are much closer to 
the known experimental data than the theoretical results. This allows to conclude that a high degree of reliability of the 
results. This is explained by the fact that the theoretical model does not take into account, for example, the length of the 
elastic shell. When the grid is thickened 2 times, the critical efforts of the models obtained with the help of the APM FEM 
KOMPAS-3D, do not practically change. Also, the forms of buckling of a cylindrical elastic shell are almost unchanged 
[ 21 ]. 


The analysis of the tables and figures allows concluding that the solutions obtained are close to the desired exact 
solution, since when the grid is thickened, the values of critical forces and forms of buckling, with the exception of some 
curvatures, for example, in 4 and 5 forms of natural oscillations of the shell, have changed [21]. 

3. CONCLUSIONS 

• The analysis of the research carried out allows the following conclusions to be drawn: 

• the Bubnov-Galarkin variational high-precision method solved in a first approximation the original nonlinear 
problem of the theory of elasticity and structural mechanics of machines by definition with a margin of bearing 
capacity of local deformations and stability of the steel shell of a roller of a road roller with varying geometry of 
the contact surface compacted until the residual displacements cease [21]; 

• developed one-dimensional physical and mathematical models of cylindrical bending, clearly illustrated with 
characteristic numerical examples of a circular shell (K) and elliptical (E) outlines [2, 21] (Figures 2); 

• in the process of the performed calculations, it was proved that the considered elastic element with a thickness of 
6.5 mm and a radius of 600 mm was deformed within the Hooke's law at stresses of not more than 86.9% of the 
yield strength tj^=2270 MPa steel grade 60C2XA [2], and theoretical safety factors =5,8 and ■?4^'=2,62 

were much higher than the minimum acceptable value [?a^r]=l,5 [8, 21]; 

• it was found that to prevent the occurrence of residual deformations, local dents and clicking [4] of the shell, it is 

JV 

necessary to use only high-quality spring-spring steel with high g >1800 MPa (- f] to make it, which is of great 
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importance in this particular case - with high requirements for fatigue strength [2, 15] and local resistance, 
corrosion and thermal resistance; in this regard, it is possible to recommend, for example, alloys on a cobalt- 
chromium-nickel base of the grades 40KXHM, 40KXHMBTU (GOST 10994-74) with = 1800 2500MPa 

JV 

(—-'1 [14-20], and also some other special structural elastically deformable steels [15] (60C2FA, 60C2XFA, 
60C2XA, 65C2BA, etc.); 

• the results of mathematical modeling and FEM analysis make it possible to optimize the physical and geometric 
characteristics of the flexible shell during its design, thereby ensuring reliable operation of the roller with an 
adjustable outline of the working surface of the drum, economical consumption of expensive spring-spring steel 
and high-quality compaction of the road surface [21]. 
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